Classical invariants and the quantum-classical link 
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The classical invariants of a Hamiltonian system are ex- 
pected to be derivable from the respective quantum spec- 
trum. In fact, semiclassical expressions relate periodic orbits 
with eigenfunctions and eigenenergies of classical chaotic sys- 
tems. Based on trace formulae, we construct smooth functions 
highly localized in the neighborhood of periodic orbits using 
only quantum information. Those functions show how clas- 
sical hyperbolic structures emerge from quantum mechanics 
in chaotic systems. Finally, we discuss the proper quantum- 
classical link. 
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For more than 20 years a lot of effort has been done in 
order to clarify the interplay between quantum and clas- 
sical behavior in chaotic systems. Periodic orbit theory 
shows the important role that classical invariants play in 
quantum mechanics. For example, the Gutz wilier trace 
formula relates the quantum spectral density with clas- 
sical periodic orbits [l]]. 

Less is known about the interplay between the neigh- 
borhood of periodic orbits (the stable and unstable man- 
ifolds) and quantum mechanics. The first step was given 
by Bogomolny |2[ who found the semiclassical contribu- 
tion of periodic orbits and their vicinities to wave func- 
tions smoothed in energy and position. Advances in this 
issue contribute, among others, to the study of the mor- 
phologies of wave functions of chaotic systems. Espe- 
cially, this gives new insights to the debated problem, 
the scarring phenomena || , an anomalous localization of 
quantum probability density along unstable periodic or- 
bits in classically chaotic systems. For example, Kaplan 
and Heller j| proposed test states constructed by coher- 
ent wave-packet sums centered on a periodic orbit. The 
constructed wave function lives not only on the periodic 
orbit, but also along the linearized invariant manifolds. 
On the other hand, Vergini and Carlo |||| developed a 
semiclassical construction of wave functions on periodic 
orbits and their neighborhoods (the stable and unstable 
manifolds) which are highly localized in energy. To re- 
semble the classical hyperbolic structure, these construc- 
tions use the orbit and the linearized dynamics around 
it. The question that naturally arises is if the classical 
hyperbolic structure is present in quantum mechanics. 
An important contribution in this direction was given by 
Nonnenmacher and Voros @. They studied the eigen- 
states structures around a hyperbolic point in a one di- 
mensional system. 



In this letter, we address that question by construct- 
ing smoothed functions living on the periodic orbit and 
their manifolds and we use only the quantum informa- 
tion (the eigenvalues and eigenfunctions) of the system 
. This clearly will show the skeleton on which the quan- 
tum mechanics is built. Moreover, some features of these 
smoothed functions confirm predictions by Bogomolny 
and are well understood with the semiclassical construc- 
tion of Ref. §. 

Our starting point is the Gutzwiller trace formula [0. 
This formula evaluates the quantal spectrum of energy of 
a bounded system in terms of the classical spectrum of 
periodic orbits. We are going to consider billiard systems 
where the classical motion is described by straight lines 
between consecutive bounces with the boundary. The 
incoming and outgoing trajectories at a bounce satisfy 
specular reflection. A periodic orbit S is characterized 
by the length Lj, and the Maslov index v§. The quan- 
tum mechanics is given by the set of eigen-wavenumbers 
k n and eigenfunctions (j) n {f) which are defined by the 
Helmholtz equation (V 2 + k 2 l )<p n (f) = , with Dirich- 
let boundary conditions. For these systems, the trace 
formula reads Hi 
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and Ms is the monodromy matrix of the orbit 5. p(k) is 
the mean level density 

The Fourier transform of Eq. ([!]) gives the classical 
spectrum in terms of the quantal one 0] 
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In particular, f(L$) — k As/2tt and \f(L)\ 2 has peaks 
at the length of the periodic orbits. Fig. |l] shows \f(L)\ 2 
obtained from the eigen-wavenumbers up to k = 264 for 
the desymmetrized stadium billiard, a classical chaotic 
system The position of the peaks agrees with the 

length of the orbits within an error (1/fc) given by the 
used window. The height of the peaks and the phases 
are also well reproduced. 
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Phase coherent contributions occurs in Eq. (||) when 
L = L$. The major contribution to the peaks is provided 
by the eigen-wavenumbers close to the Bohr-Sommerfeld 
quantization rule, k « I 25 " ( n + l) wnere n is an inte- 
ger. That is, the eigen-wavenumbers distribution is not 
random and has information of the periodic orbits of the 
system |J . The same behavior is expected for wave func- 
tions. Then, we propose the following spectral function 
for the orbit 5 

C s (f) = J2 l^(0! 2 cos(k n L s - iirvg/2). (4) 

k n <k 

We expect this function will retain the relevant features 
of the orbit contained in the eigenfunctions. In order to 
verify this point, we computed the spectral functions (Eq. 
(|])) for the orbits shown in the inset of Fig. [l] using the 
first 9910 eigenfunctions of the desymmetrized stadium 
billiard (corresponding to < k < 264). These spectral 
functions are shown in Fig. ||. A Gaussian smoothing 
with standard deviation a = 0.02 was applied in order 
to avoid the background random fluctuations. Moreover, 
with this amount of eigenfunctions, the spectral func- 
tions of all the periodic orbits up to L ~ 14 are well 
reproduced. 

These smoothed functions are clearly localized in the 
neighborhood of the orbits and some important features 
are observed. We notice that the spectral functions local- 
ize transversely to the orbits in two different ways. For 
orbits (1) and (3) of Fig. ^ the maximum of the transver- 
sal localization occurs at some distance from the orbit. 
Two peaks are observed at each side of the orbit (see 
Fig. |3| (a)). On the other hand, for orbits (2) and (4) of 
Fig. uTthe two peaks are less clear and the function £ is 
high on the orbit itself (see Fig. ^ (b)). According to Bo- 
gomolny || these features depend on the characteristics 
of the orbit but they were not specified. They come out 
clearly by comparison with the semiclassical construction 
given in Ref. (5). From that viewpoint, these functions 
are a combination of the square of the scar functions with 
even and odd transversal excitations. When the classical 
transversal motion is hyperbolic with reflection, the odd 
scar function quantize at the anti-Bohr energies. Then, 
the contribution of the odd scar functions to Eq. (||) 
change their sign, producing other mixing with the even 
scar functions (e.g. orbits (1) an (3)). 

Another important characteristic observed in the spec- 
tral functions (also predicted in Ref. ||) is their behavior 
near self-focal points of the orbit. There, the function 
is enhanced and the transversal width becomes smaller. 
An extra increase of the function was also observed at 
self-crossing points of the orbit and in the vicinity of re- 
flections with the boundary (see Fig. [2]). In a self- focal 
point the linearized stable or unstable manifold lies along 
the direction of the transversal momentum. This shows 
that the spectral functions have the imprint of the un- 
stable classical motion [frl. 



The last point should be clarified, therefore we have 
calculated a phase space representation of the spectral 
function using the Husimi distribution |ll|] of the wave 
function instead of |</>(r)| 2 m Eq.(4). As usual, we have 
used the Birkhoff coordinates. In the bottom right panel 
of Fig. [I| the phase space representation of the spec- 
tral function of the orbit (1) is showed. The stable and 
unstable manifolds of the orbit are also plotted. It is ob- 
served clearly that the spectral function lives along the 
manifolds. Furthermore, it is interesting to determine 
how this hyperbolic structure emerges in the semiclas- 
sical limit. This limit occurs when h — > or k — > oo. 
Figure || shows the phase space representation of the 
spectral function of orbit (1) using the first 150, 500, 
2000 and 6500 eigenfunctions of the system. Increasing 
the number of eigenfuntions, as we are going to the semi- 
classical limit, the functions become more localized over 
the manifolds. 

Finally, there is an important point to be considered. 
A remarkable aspect of Eq. (||) [or (||)] is that it only 
permits to obtain periodic orbits with length satisfing 
L < \n(kL)/h (with h the topological entropy). This 
is because the width of the defined peaks is ~ 1/k 
and the density of peaks increases exponentially with L 
(~ exp(ftX)/L). On the other hand, in order to recover 
the quantum spectrum with the trace formula [Eq. ([[])], 
we need all the periodic orbits with length up to the 
Heisenberg length Lh — kA (with A the area of the 
billiard). In this way, we arrive to an asymmetrical situ- 
ation where the classical information obtained from the 
quantum spectrum is wholly insufficient to recover the 
initial quantum information. This apparent contradic- 
tion was recently showed up with the development of 
a semiclassical theory of short periodic orbits ||. In 
it, the length of the required periodic orbits (to obtain 
the eigenvalues and eigenfunctions up to k) is lower than 
L < hi{L nh) / h = \n(kAh)/h. So, we can go from quan- 
tum mechanics to classical one and vice versa without loss 
of information, and the quantum-classical link is properly 
established. 

In conclusion, we have constructed highly localized 
functions in the vicinity of periodic orbits using only 
quantum information. With this evidence, we have 
shown that the classical hyperbolic structure of unsta- 
ble periodic orbits is contained in the eigenfunctions of 
the system. 
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FIG. 3. Spectral function of the orbit (3) and (2) of Fig. 1 
evaluated over the transversal section to the orbits showed in 
the insets. The transversal coordinate is y, with y = on the 
orbit. 



FIG. 1. \f(L)\ 2 evaluated from the first 9910 eigenvalues 
of the desymmetrized stadium billiard with radius r = 1 and 
straight line of length a = 1. The labels over some of the 
peaks indicate the corresponding periodic orbits shown in the 
inset. 
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FIG. 2. Spectral function (Eq. (3)) corresponding to the 
orbits (l)-(4) (from upper left to bottom right) shown in the 
inset of Fig. 1. Eq. 4 was evaluated for the first 9910 eigen- 
functions of the desymmetrized stadium billiard. 



FIG. 4. Phase space representation of the spectral function 
for orbit (1) using the first 150, 500, 2000 and 6500 eigenfunc- 
tions of the desymmetrized stadium billiard. The stable and 
unstable manifolds are plotted with points. The coordinate q 
is the arc length at the boundary normalized by the perimeter 
of the billiard and p is the fraction of tangential momentum 
P • 
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